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INTRODUCTION 

The  purpose  of  the  present  analysis  is  to  determine  the  effect  of  an  arbi- 
trary wind  distribution  upon  the  trajectory  of  the  LOKI  missile.  The  para- 
meters are  presented  in  a fashion  to  facilitate  the  design  of  a computer  to 
correct  for  wind.  Further,  the  parameters  are  adaptable  to  correcting  for 
varying  winds  which  may  exist  during  the  missile's  flight.  The  effect  of  winds 
is  treated  to  rocket  burnout  only.  The  problem  is  treated  in  two  parts;  1.  e. , 
the  calculation  of  sideways  momentum  and  angular  momentum  lnparted  by  a 
unit  wind  impulse,  and  step-by-step  calculation  of  the  resulting  flight  path  de- 
viation due  to  the  Impulse. 

A drawing  of  the  LOKI  is  shown  in  Figure  1 and  the  performance  character- 
istics are  given  in  Table  I. 

ANALYSIS 

The  LOKI  missile  rotates  at  a speed  of  approximately  100  radians  per  sec- 
ond. As  a result  there  arises  the  so-cal’ed  Magnus  effect  which  produces 
(1)  a circulation  around  the  missile  and  a resulting  lift  force  and  (2)  gyroscopic 
moments  due  to  rotation.  Both  of  these  terms  have  been  neglected  in  the  pres- 
ent analysis  on  the  basis  of  a careful  study  previously  carried  out  by  Dr.  J.  E. 
Froehlich  of  the  California  Institute  of  Technology  Jet  Fropulslon  Laboratory. 
Dr.  Froehlich  reported  that  he  had  included  these  effects  in  some  exact  LOKI 
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trajectory  calculations  involving  a constant  wind  and  found  their  contribution 
to  be  negligible.  In  the  preaent  analysis  it  ia  farther  assumed  that  the  effect 
of  thrust  malalignment  upon  the  trajectory  may  be  treated  as  a separate  phenom- 
enon, thus  may  be  neglected  here.  The  wind  is  treated  as  a perturbation 
an  the  trajectory  with  no  wind. 

The  general  equations  of  motion  in  one  plane  and  a figure  for  the  geometrl 
cal  relationships  are  as  follows: 


^AL0MG  FLIGHT  PATH 

Fj  cose*  — D -WcOS^L+X)  ' mV  - 0 

flight  path  = 0 

• °<  + Fj  3in°<  + Wsin(())i_+^)  - mV  £ = 0 

EM  =0 

16  -Mg  -6  — =0 
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The  velocity  V.  along  the  flight  path  has  bees  obtained  from  Reference  A for 
the  case  of  no  wind.  It  is  assumed  in  this  analysis  that  the  forward  speed 
will  not  be  altered  appreciably  by  winds  because  of  the  high  forward  acceler- 
ation (approximately  125  g's  Initially  and  25C  g's  at  end  of  burning)  and  high 
speed.  There  is,  therefore,  no  need  for  Including  equation  (1)  in  this  analy- 
sis and  it  has  been  dropped.  It  is  to  be  noted  that  the  third  term  of  equation 
(2)  causes  Right  path  curvature  due  to  gravity.  This  effect  is  accounted  for 
in  the  trajectory  without  wind  and  thus  may  be  neglected  in  this  wind  pertur- 
bation analysis.  A definition  of  terms  in  equations  (2)  and  (3)  is  included  in 
Appendix  A,  the  jet  damping  derivative  in  Appendix  B and  the  aerodynamic 
damping  in  Appendix  C. 

Rewriting  equations  (2)  and  (3)  using  * 8-!f  and  setting  c.o*tf  m I.O, 

3 1 n Y — and  sm=*r  — because  these  angles  are  small,  we  find 

+ Rj 


y - 


e 


m V 


(0-s) 


PH 


M 


A* 


e 


M„ 


(©-«) 


(2a) 

(3a) 


Thus  far,  the  effect  of  wind  has  not  been  taken  into  account  in  the  equation 
of  motion.  The  wind  may  be  taken  as  an  increment  in  angle  of  attack,  ^ 
where  w is  the  wind  velocity  component  (feet/second)  normal  to  the  flight 
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path  and  In  the  plane  under  consideration.  The  flight  path  angle,  '6  , during 

booster  rocket  burning  changes  less  than  two  degrees  so  the  wind  may  be  brok- 
en up  into  its  three  components  with  respect  to  the  launching  direction,  wx,  wy, 
and  wx  The  effect  of  wx  (the  component  along  the  flight  path)  is  negligible  as 
stated  previously.  The  effect  of  Wyjuid  wz  are  considered  independent  and  are 
treated  separately. 

We  now  write  equations  (2a)  and  (3a)  with  the  wind  as  a forcing  function. 


+Fj 


mV 


(e-'g) 


mV 


w 

V 


(2b) 


e - 


M 


•M 


'>e 


e-  y (e-*)=-y 


— (3b) 

V 


Integrating  equation  (2b)  with  respect  to  time,  we  get 


(4) 


Now  if  we  consider  a wind  blowing  from  time  zero  to  at  and  if  we  let  at — *0 
adille  holding  the  product  w-at  * 1.0  so  that  we  have  a unit  wind  impulse,  we 
get,  upon  completing  the  integration, 


L< 


A*i(t>  " “ tv  V 


(5) 
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0 

aY  i (t)  then  represents  the  immediate  effect  of  the  wind  impulse  upon  If. 

Similarly  we  obtain 

*Qi  (t)  = — • Y (6) 

from  equation  (3b). 


and  . -L.  are  considered  constant  in  the  above  integrations  be- 
cause At * 0.  It  is  to  be  noted  that  A%i(t)  is  independent  of  forward 

velocity.  This  is  shown  by  substituting 


L * > V 1 b 1 
L°<  2 i'  d <x 


in  equation  (5) 


t) 


'zips/1  b1  f 


* dC^ 

<X 


m 


1/  _ U * 


rn 


(7a) 


A©  i (O  increases  as  a function  of  V and  is  shown  by  rubstitutlng 
*<*-  '/iPV'b5C^  in  equation  (6) 


Aewt)  = ^^b3c- 


i 


v 


(7b) 


Curves  of  afyft)  and  nG(  (t)  are  shown  in  Figures  2 and  3.  a^(t)  and 

i 

A©(,  (t)  are  then  the  increments  developed  immediately  upon  application  of 
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the  unit  wind  impulse.  To  determine  the  error  at  the  end  of  rocket  burning,  an 
Integration  must  be  made  recognizing  the  fact  that  the  coefficients  vary  with  time. 
Equations  (2a)  and  (3a)  were  rewritten  in  the  form 

y * 2 a,  (9-K)  » 


where 


0 = - 

+ c0o 

z>,  * 

La*  ■+•  Fj 
nr\  V 

Z) i.  - 

[Md4 

I 

w „ 

- 

1 

Using  c*<  * 0 - , the  simultaneous  solution  of  equations  (8)  and  (9)  becomes 

*^64  -i-2. ^ i + co,,* ) •=".  — C 

^ e = , *►  ^ o 


then  letting 
and 


^>e  = *4”  ?!«  ^ 2.^  , •+ C00 


equation  (10)  becomes 

C ik  2.  e -f0O^cK=O 


As  a solution  let 


= Q,  £ 4i  cos  cot  -tB  sm  uA}’ 


where 


CO  = yco**  “ 
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From  equation  (12)  we  obtain 

.-V 


6.  Ae^A cot  wt  ■♦’B*'*'  u>t}  •+•  {-Ago  -♦-Bojcos cot^J 


(13) 


From  equatlona  (12)  and  (13)  we  find  att«o 


°<o  ~ * 

•< « — Bco  — )ie  A “ Bw  — ?\e«x>0 
B £ 

Thus  equatlona  (12  and  (13  become 

6 * ' ^«oCOS<*>t,  «♦*  L~~?)  5**  u>t,\ 

«*  (t.)  — + Slh  u>t,-K*#Co*u*t,] 

From  equation  (8)  t t 

yet,)  ~ +J  z^.ch at  * *^-*-5^  • At 

A sample  step -by -step  calculation  using  the  above  method  was  run  for  a 


(14) 

(15) 

(16) 


gust  occurring  at  . 3 of  a second  after  Initiation  of  burning.  This  was  run  as 
a check  point  on  the  work  done  by  Dr.  Froehlich  and  is  shown  in  Appendix  D. 
This  check  agreed  with  Dr.  Froehllch's  calculations  within  five  percent. 


An  attempt  to  compute  the  trajectory  step-by-step  using  a linear  extrapo- 
lation over  the  time  interval  proved  to  be  unsatisfactory,  due  to  the  very  high 
accelerations.  The  above  described  sinusoidal  extrapolation  was  more  labor- 
ious but  gave  satisfactory  accuracy  with  manual  computation. 

i 

of  2 2\  , and6j,  are  given  in  Figs.  4,  5,  and  6 and  in  Table  II. 
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The  final  values  of  flight  path  deviation  at  fuel  burnout  <*X  f (i)  due  to 
the  wind  impulses  are  shown  in  Fig.  7 plotted  versus  time.  (The  values  for 
this  curve  were  determined  from  work  done  by  Dr.  Froehlich. ) This  term 
aXp  (t)  which  we  shall  call  the  influence  coefficient,  is  then  the  integrated 
effect  of  a unit  crosswind  impulse  applied  at  a particular  time. 

The  total  flight  path  angular  deviation,  aXt  , due  to  a continuous  dis- 
tribution of  winds  may  be  calculated  by  integrating  the  product  of  the  influence 
coefficient  and  the  crosswind  in  feet  per  second  from  zero  to  the  time  at  the 
end  of  the  boost,  the  equation  being 

= / [oXf  (t)  ' w(t)]  dt  (17) 

where  aXf  (-t)and  w<ft)  are  the  values  of  the  influence  coefficient  and  the 
wind  at  the  time  under  consideration. 

It  is  to  be  remembered  that  the  value  of  X given  by  equation  (16)  is  the 
direction  of  the  flight  path  of  the  center  of  gravity  of  the  complete  LOK3.  The 
direction  of  the  c.g.  of  the  DART  (the  forward  body)  differs  from  this  because 
of  the  distance  between  these  two  c.g.  's  if  there  is  a residual  0 . X p does 
not  completely  define  the  eventual  direction  of  the  DART,  and  X are  also 
boundary  conditions.  Because  of  lack  of  data,  it  is  not  possible  to  evaluate 
these  two  additional  effects,  however,  when  data  becomes  available  an  evalua- 
tloaof  them  will  be  made.  They  will  be  most  effective,  of  course,  wheh  gusts 
occur  late  in  the  boost  phase. 
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Tc  facilitate  the  determination  of  the  error  due  to  wind  when  the  wind  is 
known  as  a function  of  distance  along  the  flight  path  instead  of  in  terms  of 
time,  we  shall  use  a wind  impulse  defined  as  W • — 1.0  . On  this  basis 

we  may  show  equations  (2b)  and  (3b)  as 

/**•? l/N#® '*>"•&  /-V7 


©dt 


ds 

d t / I 


**-&-[ T(0-Wl1f-/ T V'ds 


however  — * V,  so 
d t 


>*  -jm  °‘*  S' T »-"*  /t'T7  * 

Now  integrating  as  we  did  to  get  equations  (5)  and  (9)  and  substituting  for 
w-as  “ 1.0  , we  get 


^ ■ V* 


(18) 


M-  I 

&©t(S)—  j.  ' ■^~3L 


(19) 


It  is  seen  that  ^ and  o©i,(s)  — A0<(t)  so  that 

4.Y^(s)  decreases  as  a function  of  V and  & ©</*)  is  independent  of  V. 


CONFIDENTIAL 
Security  Information 


CONFIDENTIAL 
Security  Information 

-10- 

Corves  of  and  are  shown  as  Figs.  (8)  and  (9).  The  sharp  peak 

(negative)  in  &&{,($)  at  180  feet  is  due  mainly  to  the  relatively  high  static 
stability  margin  at  this  point.  The  static  margin  is  approximately  three  times 
as  high  at  this  point  as  at  other  points  in  the  trajectory.  The  rise  in  o0i,(S) 
towards  the  end  of  burning  is  due  to  the  decreasing  moment  of  inertia.  Fig.  10 
shows  the  influence  coefficient,  ( S)  , as  a function  of  distance  where 

(S)  = ^-  ■ a*F<t). 

The  total  flight  path  angular  deviation  due  to  winds  may  thus  be  obtained  by 

s 

=J  (S>-W(s)]  ds  (20) 

o 

or  by  use  of  equation  (18). 
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TABLE  I 

LOE  PERFORMANCE  CHARACTERISTICS 
Initial  weight  of  booster  plus  Dart  « 24.  30  lbe. 
Final  weight  of  booster  pins  Dart  * 12. 23  lbs. 
Weight  of  booster  propellant  ■ 12. 07  lbs. 

Weight  of  Dart  only  ■ 8. 3 lbs. 


! CHARACTERISTICS  AS  A FUNCTION  OF  TIME 

1 

Time, 

Thrust, 

1CG 

Mass 

Velocity 

Seconds 

Pounds 

m = i/g 

Ft/sec 

0 

0 

4.14 

.756 

0 

.10 

3120 

4.00 

.718 

365 

.20 

3200 

3.84 

.673 

825 

.30 

3275 

3.65 

.627 

1300 

.40 

3300 

3.45 

.580 

1830 

.60 

3326 

3. 24 

.534 

2400 

.60 

3300 

2.99 

.487 

3025 

.70 

3160 

2.74 

.440 

3690 

.80 

2125 

2.49 

.400 

4260 

.90 

620 

2.36 

.382 

4475 
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STABILITY  AND  FLIGHT  PATH  CHARACTERISTICS 


Time, 

Seconds 

USA 

egg 

mm 

mm 

-1. 16 

0 

.0 

-1.15 

.575 

.0 

-1.20 

- .0916 

- 1.056 

-1.292 

27.519 

.646 

1.030 

-1.25 

- .272 

- 6.76 

-1.  522 

12.08 

.761 

2.60 

-1.  275 

- .437 

- 18.4 

-1. 712 

8.  01 

.856 

4.290 

-1.  315 

- .650 

- 40.8 

-1.  965 

6. 19 

.9825 

6.36 

-1.  365 

- .950 

- 108.0 

-2.315 

5.30 

1.157 

10.40 

.30 

-1.44 

-1.28 

- 234.0 

-2.72 

4.84 

1.330 

15.28 

.35 

-1.  60 

-1.45 

- 302.0 

-2.95 

4.  51 

1.475 

17.40 

.40 

-1.60 

-1.71 

- 348.0 

-3.31 

4.28 

1.655 

18.66 

.46 

-1.68 

-2.137 

- 404.0 

-3.817 

4.22 

1.908 

20. 16 

.50 

-1.80 

-2.50 

- 441.6 

-4.30 

4.25 

2.150 

21.04 

.55 

-1.  92 

-2.89 

- 480.4 

-4  81 

4.37 

2.405 

21.96 

.60 

-2.07 

-3.34 

- 553.2 

-5.41 

4.48 

2.705 

23.52 

.66 

-2.23 

-3.79 

- 671.2 

-6.02 

4.81 

3.  01 

25.  96 

.70 

-2.40 

-4.25 

- 802.8 

-6.65 

4.89 

3.325 

28.32 

.76 

-2.60 

-4.74 

-1020.4 

-7.34 

5.  01 

3.67 

31.98 

.80 

-2.84 

-5.27 

-1234.0 

-8.11 

4. 83 

4.  055 

35.10 

.86 

-3.  00 

-5.82 

-1474. 0 

-8.82 

4.49 

4.410 

38.40 

.90 

-3.08 

-6.40 

-1692. 0 

-9.48 

4.35 

4.74 

41.10 
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TABLE  m 


INFLUENCE  COEFFICIENTS  (t)  and 

a*,.  (5)  = 


w*t«  1.0 
at  Time, 
sec. 

Radians 

V 

Ft/Sec 

Radians 

2> 

Diet-,  Ft. 

0 

0 

0 

0 

.05 

-0.56 

150 

-3.70 

5 

.1 

-1.04 

365 

-2.85 

15 

.2 

-1.05 

825 

-1.27 

0 

74 

.3 

-2.31 

1300 

-1.78 

179 

.4 

-1.30 

1830 

-0.71 

337 

.5 

-1.01 

2400 

-0.42 

547 

.6 

-0.91 

3025 

-0.  30 

818 

.7 

-0.18 

3690 

-0.05  . 

1154 

.8 

-0.64 

4260 

-0. 15 

1553 

.9 

+0.76 

4475 

+ 0.17 

1995 
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APPKNDIX  A 
LIST  OF  8YMBOLS 


Lm 


W 

y 

m 

t 

I 


dL 


\ 

T 


p V • b 


d Cl 

d o< 


, slope  of  the  lift  curve  in  pounds  per 
radian 


p-  . 00230  aix-  density  In  slugs/cu/ft. 

V > forward  velocity,  ft/sec 

b a diameter  of  booster  body,  (b  a . 25  ft) 

dCi.  J /" 

= slope  of  lift  coefficient  curve  (assumed  equal  to  ~ from 
dot  a °\ 

Reference  (A),  the  slope  of  the  normal  force  coefficient  curve) 

per  radian 


o<  m angle  of  attack  of  missile  body  axis,  radians 
a rocket  jet  thrust,  pounds  (obtained  from  Reference  A) 

* instantaneous  weight  of  missile,  pounds 

« angle  between  vertical  and  tangent  to  flight  path,  radians 

* instantaneous  mass  of  the  missile,  slugs  (obtained  from  Reference  A) 
q.  radians  per  second 

. instantaneous  moment  of  inertia  of  the  missile  about  the  Y axis, 
slug -ft?  |m  x(K^)2  . I<4k,  obtained  from  Reference  A). 


* d£ 
“ I c.o. 


i 
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APPENDIX  A (contd) 
LIST  OF  SYMBOLS 


d*0 

dt* 


* radians  /sec2 


fS  - ^ + It  ~ + M»i 

rotational  damping  due  to  rocket  jet.  The  derivation  of  this  term  is 


shown  in  Appendix  C . 


"a*  - rotational  aerodynamic  damping.  Derivation  shown  in  Appendix  C - 

k/t  6 M 1 w*  I* 

* TS  =T  Pv  b Cm, 

\ c&  » distance  from  aft  end  of  booster  to  Instantaneous  center  of  gravity, 
calibers  (obtained  from  Reference  A) 

\cp  - distance  from  aft  end  of  booster  to  Instantaneous  center  of  pressure, 
calibers  (d  xlp  from  Reference  A) 

C-v-dc-lc*)^ 

e - 
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APPENDDv  B 
JET  DAMPING  MOMENT 

The  effect  of  rocket  jet  damping  due  to  angular  missile  velocity  was  ob- 
tained from  Reference  (b).  The  jet  damping  term  of  equation  (3)  from  this 
reference  is: 


r v*  M j . 


e 


where  Q =0  = missile  airframe  angular  velocity  - radians/sec 
incremental  mass  of  gases  at  any  section 
distance  of  particle  of  gas  under  consideration  from  eg  - ft. 
velocity  (ft/sec)  of  gas  particle  under  consideration  with  respect 
to  missile 


r = 

V = 


Now  pj  Av  * c,  where 

C - mass  flow  in  slugs  per  sec.  at  a particular  station 

A - cross  sectional  area  at  section  under  consideration  ft2 


mass  density, 

c 

"Hr 


er 

V = 


sec* 


fi 


ft 


M 


j,  ^Z^irvypjAdc-r-p  -z/erdr 


CONFIDENTIAL 
Security  Information 


CONFIDENTIAL 
Security  Information 

-17- 


APPENDIX  B (Contd) 

JET  DAMPINQ  MOMENT 

Assuming  equal  rate  of  burning  all  along  booster  and  neglecting  nozzle  length 
a plot  of  "c"  would  be: 


r 

and  jCdr  — 5 * • = 1.15.  The  integral  of  Crdr  thus  equals  1.15  l. 

J 2 

The  distance  of  the  centroid  olj Cd  r from  the  tail  end  of  the  booster  » 
1/3x6  ft  * 1.867  ft. 
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AP PEND IX  B (Contd) 
JET  DAMPING  MOMENT 


<11 

(2) 

<31 

<41 

(5) 

(6) 

<71 

Diet,  from  CO  to 

Zrr\i  r v 

gegen 

K* 

Time 

Tail  End  of  Booster, 

P4tt 

* fc.t6t 

1 

I CO 

feet 

• (*)-!. 647 

= 1.1*0) 

sn 

I 

3.  75 

2.  08 

2.39 

4.  78 

4. 14 

1.  154 

3.80 

2. 13 

2.45 

4.90 

4.00 

1.25 

3.87 

2.20 

2.53 

5.  06 

3.84 

1.319 

3.97 

2.30 

2.64 

5.28 

3.65 

1.445 

4.07 

2.40 

2.76 

5.52 

3.45 

1.60 

.5 

4.20 

2.53 

2.91 

5.82 

2.24 

1.798 

.0 

4.35 

2.68 

3.08 

6. 16 

2.99 

2.06 

pas 

4.54 

2.87 

3.30 

6.60 

2.74 

2.41 

WSm 

4.73 

3.06 

3.52 

7.04 

2.49 

2.83 

4.83 

3. 16 

3.63 

7.26 

2.36 

3.08 
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APPENDDC  C 

AERODYNAMIC  DAMPING  IN  PITCH 


The  damping  moment  derivative  which  is  denoted  by  ■ is  obtained  in 

the  following  maimer: 

Consider  the  missile  to  have  a lifting  surface  i located  at  a distance  x; 
from  the  c.  g. 



V 


This  surface  will  experience  an  induced  angle  of  attack  eq’.v»l  to 
to  the  angular  velocity  0 in  pitch.  The  Moment  M*  J ^ ' 
acts  in  a direction  opposite  to  that  of  0 . The  total  moment  is 
summing  over  all  the  lifting  surfaces. 


The  damping  moment  derivative  which  is  obtained  by  differentiating  the  above 
expression  with  respect  to  6 is 


xi  8 
V 

V 

obtained  by 
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